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8.2 MISREGISTRATION'S EFFECTS ON CLASSIFICATION *D PROPORTION ESTIMATION 
■r ACCURACY* ! 

R^. Juday and F. Hall, Johnson Space Center 


The estimates of crop type and acreage (proportion estimation) are undertaken 
in the AgRISTARS program by registering multiple date acquisitions of small 
subareas of Landsat scenes (termed segments), and applying mul tispectral anal- 
ysis to them. An important contribution to errors in classification and acre- 
age estimates is misregistration between multiple acquisitions. The relation- 
ship can be expressed as: 
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The particular series of operations applied are shown diagrammatical ly in 
Figure 1. The spectral data brought in with ancillary data are transferred 
into the Kauth-Thomas space, where features are extracted. The Kauth-Thomas 
greeness is a function of time for multiple acquisitions. It is approximated 
in functional form with various parameters and the vector of these parameters 
and its distribution within the segment you are working with is decomposed so 
that its overall distribution function is approximated a posteriori . The a 
priori probabilities given class and spectral vector that describe a class are 
the parameters that give you the estimate of the area and crop type. The an- 
cillary data is the place where you bring in the fact that you are working in 
a specific area, thereby limiting the variety of possible crop choices. 


Figure 2 illustrates the taking of a Landsat feature vector and deriving the 
brightness and the greeness. A procedure derives certain pa-aneters that de- 
scribe the time at which a crop begins to green up, and the steepness with 
respect to which the curve moves off the soil line is a function of time; time 

being a characteristic difference between some of the crops that we're trying 
to differentiate. This kind of procedure is sensitive to absolute radiometry 
to perhaps a larger extent than some researeners maybe are interested in with 
regard to maintaining radiometric fidelity in the resampling process. Figure 
3 illustrates the procedure whereby a scene is broken into various strata. 
There are distribution functions within each strata for the spectral measure- 
ment vectors via a procedure called CLASSY. CLASSY breaks the curve down to 
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several distributions which are assumed unimodal and normal. This is the 
place at which the proportion estimates fall out. In the process of doing 
this, the mixture pixels along the boundaries of the segments create confu- 
sion. Registration error looks something like a mixture pixel. If there is 
misregistration, and a given address for a pixel jumps across the boundary. It 
acts as a noise, and in the multitemporal sense will tend to confuse crop dis- 
crimination. 

Figure 4 presents the expression of pure and mixed pixel contributions to 
classification error. Knowing full wll that we have mixture pixels, and if 
you can parameterize the problem in terms of your estimate of a crop that has 
found its truer pixels, the equation can be applied to derive a final propor- 
tion estimate. The variance in the estimate of the crop proportion in mixed 
pixels is typically larger than the variance cf Ihe estimate of the crop that 
is in pure pixels. You would like to have q, which is the proportion of the 
scene which is in pure pixels, to be a larger number so as to reduce the over- 
all variance. That, of course, is done in one way by going to smaller pixels 
and in another by having better registration accuracy so that you don't have 
the multitemporal jumping of pixels. Now the quantity q will approach a num- 
ber that is limited by the spatial resolution of the system as the registra- 
tion error declines. There is a point at which you cannot drive q to be a 
larger number with increasing registration accuracy and that declining gain 
determines the point at which for any particular analysis technique, it is not 
practical to expend more effort in trying to get better. 

For any given sensor IFOV geometry, you can derive ty; al populations of 
fields and you can plot histograms of the number of fielc against field size 
and you can draw different histograms for the different K nds of crops that, 
according to your ground truth, go into making up the population examined. 
Dave Pitts and Gothem Badhwar have published some results for representative 
samples in the growing regions of the U.S. that are of significance to us 
(Figure 5). As a function of the resolution element, the IFOV of the sensor, 
according to those kinds of histograms in Figure 5, can draw the proportion of 
pure pixels in a given crop (Figure 6). This is the quantity which you would 
like to drive up. The TM will have a resolution that is considerably smaller, 
so the proportion of pixels that are pure in any given area will be larger. 
You can draw an analogy between the misregistration error as acting much the 
same as an increased resolution IFOV and expect that things will go in the 
same direction although we have some doubts that you can draw an exact equi- 
valence there because there is some difference in behavior. 

This is a workshop and not a symposium so I am not presenting a final result. 
At this point we do have a program underway to come up with quantitative num- 
bers for the results for the effect of the misregistration error. Also, we 
are tending to work mostly on the multitemporal aspects of misregistration at 
this point. Band-to-band misregistration has similar effect, but, in our par- 
ticular application, it's a less dynamic kind of an affect. We built a regis- 
tration system af JSC to a point when it was clear that the one pixel regis- 
tration accuracy that the LACIE processor was giving us was only marginally 
adequate for the LACIE technology, and with AgRISTARS we needed something 
better. The MDP, being a full-frame system, does not concentrate as heavily 
its computation power on an individual pixel. We set ou selves a goal of 2/10 
pixels as the registration accuracy that we would like to get to. This was 
based on the arbitrary criteria that what you would like is 50X of the multi- 
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temporal energy to come from a common area on the ground. When you do that, 
it’s really amazing how strictly that ties down the registration accuracy that 
you think you would like to get to. When you begin talking about registration 
at that accuracy you uncover all sorts of things that you didn’t have to think 
about previously. As Figure 7 illustrates, the difficult problem to determine 
is the registration accuracy that you want analytically when you overlay one 
aperture over another. For pixel A and pixel B. the portion of the clear 
aperture is a function of the position of the pixel B, given pixel A is fairly 
easy to calculate. But when you begin to add other pixels that complicate? 
your clear aperture upon adding a third aperture, given a position of the 
first two is some value on the ramp. It is highly nonlinear and this gets 
worse when you go to two dimensions, and when you get nonsquare IFOVs, and 
when you have a variation of response within a IFOY. Our plan is to pursue a 
Monte Carlo approach in assuming a misregistration error and model its distri- 
bution through the process. 

In conclusion, it should be noted that we do not have a quantification of mis- 
registration's effect on our analysis at levels better than we currently 
achieve. Synthetic data is one approach; extrapolations obtained by deregis- 
tration is another. Registration accuracy requirements are difficult to spe- 
cify a priori . ”Nore is better"; everyone’s on expe»*t. Extraoolation to the 
performance expected from better registration is risky, and highly dependent 
on analysis techniques. The proper arena for testing is one's own applica- 
tions analysis. 
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Figure 1 . Crop ID/Proportion Est. 
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Figure 2, Data Transformation/Fcature Extraction 
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OiGiTAL IMAGE SCENE 
Figure 3. Mixture Decomposition Estimation 
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P = qPo + (l-q)Pm 

P = PROPORTION ESTIMATE FOR SCENE FOR CROP 
Pq = ESTIMATE OF CROP PROPORTION IN PURE PIXEL AREA 
Pm = ESTIMATE OF CROP PROPORTION IN MIXED PIXEL AREA 
q = PROPORTION CF SCENE IN PURE PIXELS 

VAR P = q2no+ (l-q)2om2 + 2q(l-q)aom 

= VAR (Pq) 

“ VAR (Pfn) 

Oom “ COV (PoPm) 

Figure 4. Pure and Mixed Pixel Error Contributions 
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FICUHE - Field si/.c disUibulion for vvinlcr wheat, coni, soybeam. 
fpnruj wheat, winter wl^cat indutimg strip fallow, and spri’ig wncat 
Including strip fallow. 



CONVOLUTION OF TWO APERTURES: 
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